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ABSTRACT 

We investigate the dependence of the complete system of 22 Lick indices on 
overall metallicity scaled from solar abundances, [^], from the solar value, 0.0, 
down to the extremely-metal-poor (XMP) value of -6.0, for late-type giant stars 
(MK luminosity class III, \ogg = 2.0) of MK spectral class late-K to late-F 
(3750 < Teff < 6500 K) of the type that are detected as “fossils” of early galaxy 
formation in the Galactic halo and in extra-galactic structures. Our investigation 
is based on synthetic index values, I, derived from atmospheric models and syn¬ 
thetic spectra computed with PHOENIX in LTE and Non-LTE (NLTE), where 
the synthetic spectra have been convolved to the spectral resolution, R, of both 
IDS and SDSS (and LAMOST) spectroscopy. We identify nine indices, that 
we designate “Lick-XMP”, that remain both detectable and significantly [^]- 
dependent down to [^] values of at least ~ —5.0, and down to [^] ~ —6.0 in hve 
cases, while also remaining well-behaved (single-valued as a function of [^] and 
positive in linear units). For these nine, we study the dependence of I on NLTE 
effects, and on spectral resolution. For our LTE / values for spectra of SDSS 
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resolution, we present the fitted polynomial coefficients, Cn, from multi-variate 
linear regression for I with terms up to third order in the independent variable 
p airs (TpfF, iff]), a nd (V — K, [^]), and compare them to the htted values 


of lWorthev et al\ (119941 ) at IDS spectral resolution. For this htted / data-set we 
present tables of LTE partial derivatives, ^|[m/h], ^[m/h] l^eff, apyf^l[M/H], and 
\(v-K}, that can be used to infer the relation between a given difference. 


dl 


a[M/H] 


AJ, and a difference ATefj or A(I/ — K), or a difference A[^], while the other 
parameters are held hxed. For Fe-dominated Lick indices, the effect of NLTE is 
to generally weaken the value of / at any give Teg and [^] values. As an example 
of the impact on stellar parameter estimation, for late-type giants of inferred 
Teff > 4200 K, an Fe-dominated / value computed in LTE that is too strong 
might be compensated for by inferring a Teg value that is too large. 


Subject headings: stars: atmospheres, fundamental parameters, late-type 


1. Introduction 


The determination of stellar parameters, especially overall metallicity (denoted here 
[^] unless otherwise indicated) and detailed abundances of individual metals, in stars of 
remote Galactic and extra-galactic structures has become crucial to the study of galaxy for¬ 
mation and evolution, including that of the Milky Way. Metallicity distribution functions 
for galaxies and globular star clusters (GCs) reveal information about multiple populations 
from multiple star formation episodes and allow the investigation of the history of star 
formation. Ghemical tagging of stellar populations allows the investigation of the link be¬ 
tween Galactic structures such as GGs and nearby extra-galactic ones such as ultra-faint 
dwarf (UFD) and dwarf spheroidal (dSphe) satellite galaxies, a nd the process through w hich 
galaxies are assemble d in hierarchical s t ructu re fo rmation (see Mucciarelli et al. (120131) for 
a recent example, and Frebel fc Norris ( 20131) and Belokurov (2013) for reviews). However, 
stars in remote structures are often only signihcantly and efficiently detectable with low- 
to moderate-resolution spectroscopy, such as that of the Sloan Digital Sky Survey - Sloan 
Extension for Galactic Understanding and Exploration (SDSS-SEGUE), that precludes the 
measurement of individual spectral lines, and usefully accurate [^] values must be obtained 
from obser v ationa ll y expensive follow-up spectro scopy at high spectral resolution {e.g. see 
Aoki et al\ (120131) ). Ilvezic. Beers &: Juric (120121) contains a recent authoritative review of 
how modest resolution spectroscopic surveys are revolutionizing our study of Galactic popu¬ 
lations and leading to insights into Galactic formation. As a result, there is active interest in 
novel methods for extracting [^] and values from low- to moderate-resolution data (see 




































3 


Schlaufman &: Casevi (120141) for a recent investigation of a method based on IR molecnlar 


bands, and iMiller (120151) for a very recent investigation of a photometric method based on 
SDSS ugriz photometry.) 


The original system of 11 (IGorgas et fli.lll993l ). and then 21 (IWorthev et all (Il994l ). W94 
henceforth), Lick/IDS spectral indices, /, was defined to optimize the determination of [^] 
and, hence, age (t), from integrated light (IL) spectra of faint spatially nnresolved old stel¬ 
lar populations (—1 < [^] < -1-0.5, 0.5 < t < 20 Gyr) dominated by G and K stars, 
obtained at low spect r a reso lution (R < 1000) in the A4000 to 6200 A spectral band. 
fjWorthey fc Ottavianil fjl997l) explored four new indices based on two dehnitions, each, of 
Hy- and H(5-centered features, but found them to be of limited usefulness - a result con¬ 
sistent with our own investigation.) Because of the emphasis on old stellar populations, 
such as that of globular clusters (GCs) and “early-type” galaxies, red giant (RGB) stars 
are important contributors to the IL spectrum. The indices were discovered by empirically 
identifying composite spectral features in low R spectra of Galactic and GC G and K stars 
that showed a significant and useful correlation with one of [^], Teg, or log^f while (hope¬ 
fully) depending less sensitively on the other two. The original Lick/IDS system was defined 
with spectra obtained with the Lick Observatory Image Dissector Scanner (IDS) having a 
spectral sampling, AA, of ~ 8 A in a region centered at 5200 A spanning a AA range of 
2400 A, corresponding to spectral resolution R = A/AA of ~ 650. 


Given the usefulness of the Lick indices fo r modern mode r ate re solution spectroscopic 
surveys such as SDSS-SEGUE and LAMOST, Franchini et al. ( 2010 ) developed the system 
further by creating a synthetic library of I values for dwarf and giant stars derived from 
synthetic spectra that had been convolved to the higher SDSS R value of 1800, and tested 
the predicted relationship between the I values and stellar parameters against that of several 
empirical spectral libraries, including the SDSS-DR7 spectroscopic database itself. They also 
supplemented the 21 indices of W94 with a new near-UV index, namely CaHK - a prominent 
spectral feature that has proved us eful for identifying candidat e XMP stars in productive 
surveys such as the HK surv e y of Beers. Preston fc ShectmanI (119921) . One of the main 


conclusions of iFranchini et al\ (120101 ) is that Tefj values derived from fitting their synthetic 


indices to SDSS-SEGUE spectra of late-type giants were systematically lower than the Tgff 
values derived with the SEGUE Stellar Parameter Pipeline (SSPP). 


Many of the Lick indices are dominated by, or have signihcant contributions from, lines 
of Fe I. Initial investigations of metal-poor stars in which the Fe I extinction is treated 
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with Non-local Thermodynamic Equilibrium (Non-LTE, NLTE) indicate that NLTE effects 
become increasingly important with decreasing metallicity, and alter the inferred Tefj values 
by as much as 400 K and the derived abundai ices of Fe by 0 . 4 dex for metal-poor gian ts 
with par ameters based on R AVE survey spectra fISerenelli et a/.l (120131) . iRuchti et all (120131) 1. 
Recently IShort et al\ (1201311 and papers in that series have found that when most of the light 
and Fe-group metals that contribute to the visible band spectral line blanketing of mildly 
metal-poor RGB stars is treated in NLTE, the value inferred from spectrophotometric 
spectral energy distributions (SEDs) of R ~ 100 is reduced by > 50 K. Therefore, calibration 
of the Lick /([^]) relation based on NLTE modeling may be crucial to using the indices for 
accurate [^] inference. 


Our goal is to extend an analysis of the detectability, and sensitivity to stellar parameters, 
including [^], of the Lick indices to the regime of extremely metal-poor (XMP) red giants, 
and to investigate the magnitude of NLTE effects on the value of modeled Lick indices. In 
Section [2] we describe the atmospheric models and the spectrum synthesis, and the procedure 
for producing synthetic I values; in Section [3] we identify the Lick indices that remain most 
useful at XMP metallicities, and provide useful polynomial fits and partial derivatives for 
index values, /, modeled in LTE, in terms of Tgg, V — K, and [^]; in Section H] we present 
conclusions. 


2. Modeling 
2.1. Model grid 


We have used PHOENIX V. 15 (iHauschildt et a/.l 119991) to computed a grid of atmo¬ 
spheric r nodels and corresponding synthetic spectra in both LTE and “massive multi-species 
NLTE” ( Short fc Hauschildt 2009 ) for very- and extremely-metal-poor (VMP and XMP) red 
(and “orange”) giant stars of MK spectral class late F to late K covering the range of T^g 
and [^] value of stars that are spectroscopically accessible a t Galactic halo dis tances and 
that serve as useful stellar “fossils” for Galactic archeology (IGohen et a/.l 120131) . The grid 
also includes red giants of higher [^] value representative of the solar neighborhood and disk 
population for comparison. The parameters of the LTE grid are: 


3750 < Teff < 6500 K with AT^g = 250 K, 

-6.0 < [f ] < 0.0 with A[f ] = 1.0 for [f ] < -2.0 and 0.5 for [|] > -2.0, 
log^ = 2.0 
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Fig. [T] shows the TKin(Ti 2 ooo) structure of a subset of our models for Teg = 4000 K and 
[^] values of -0.5, -2.0, -4.0, and -6.0, where ri 2 ooo is the monochromatic continuum optical 
depth at 12000 A and serves as our standard radial depth variable. The reduction in 
the well-understood back-warming and surface cooling effects c aused by line ex t inctio n as 
[^] decreases is readily noticeable. For comparison, the grid of iFranchini et al\ (120101) has 
3500 < Teff < 7000 K with ATg^ = 250 K, —2.5 < [A] < 0.5 with A[^] = 0.5 generally, 
with the addition of [A] = —4.0 for their a-enhanced models, and 0.5 < \ogg < 5.0 with 
A log (7 = 0.5, where [^] denotes the scaled abundance parameter for elements other than 
a-process elements. The most important distinguishing features of our grid are the extension 
to [^] values of -6.0, and the inclusion of NLTE models for a subset of parameter values 
spanning the grid. Although there are few stars of [^] < —3.5, even among halo stars useful 
for Galactic archeology, extending the grid to [^] = —6.0 allows us to anchor the /([^]) fit 
through the useful [^] range. 


Because of the large number of [^] values (nine) and doubling of most of the grid to 
include NLTE counterparts, the number of models was limited by fixing the log^f value 
at 2 . 0 , representative of the giant population, and giving all models a scaled-solar abun¬ 
dance distribution based on the abundances of iGrevesse fc San vail (1199811 (therefore, for our 
grid the over a ll met allicity parameter, [A]^ is identical with the [^] parameter). W94 and 
Gorgas et al\ (Il993l) found that important Tgfj- and [^]-sensitive I values were least sensi¬ 


tive to log 77 . Nevertheless the lack of the log 77 dimension in our I polynomial fits described 
below leads to polynomial fitting coefficients that are not directly comparable to those of 
W94 from their analysis of observed IDS spectra. Adopting an a-enhancement of 0.0 is 
expected to have a minor effect on the differential comparison of most Fe-dominated [^]- 
sensitive / values computed in LTE and NLTE (because these are really more [^] indicators, 
and Fe is not an a-element), and will allow direct assessment of pure NLTE effects on I val¬ 
ues across the entire range of [^] value. Moreover, incorporating a-enhancement requires a 
model for how the value of the enhancement increases with decreasing [^] value in the range 
0.0 to -1.0. Nevertheless, a future direction is to extend our grid in the a-enhancement di¬ 
mension. As a preliminary assessment of the effect of a-enhancement, in Tabled] we present 
the computed values of the index, /, computed at SDSS spectral resolution for our subset 
of nine “Lick-XMP” indices (see below) for LTE models of Teg = 4000 K, log 77 = 2.0, for 
which the I values are relatively strong, and select [A] values of -2.0 and -4.0 and scaled 
solar abundance, and abundances with the maximal relative enhancement of -1-0.4 of the 
eight a-process elements of even atomic number from O (Z = 8 ) to Ti {Z = 22). Mg is 
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an a-process element, and the Mgi and Mg b indices are more strongly affected, increasing 
by a factor of ~ 1.5 at [^] = —2 where their value is still large. Of our nine Lick-XMP 
indices, the results presented below for Mgi and Mg b should be regarded as most suspect 
and require follow up investigation with a full a-enhanced NLTE and LTE model grids. 
The inclusion of NLTE effects in the modeling is not likely to change the conclusion that 
non-a-element spectral features are significantly less affected by a-enhancement than are 
a-element spectral features. As might be expected, the effect of a enhancement on the Fe- 
dominated indices Fe4531, Fe5015, Fe5270, Fe5335, and the Na D index is minor, being of 
the order of 10% or less. We note that for some non-o-element indices, the effect 
of o-enhancement can be to reduce the index slightly at some [^]-values, while 
increasing it at others. Changes to the composition can affect the spectrum in 
the bracketing pseudo-continuum windows that define any index, and will also 
indirectly affect the strength of features contributing to any index through the 
effect on the electron number density. 


Our models have spherical geometry with radii based on an adopted mass of M = IMq, a 
microturbulence broadening parameter of = 2.0 km s“^, which is consistent with what has 
been measured and adopted for late-type giants generally, and a mixing-length parameter 
for the treatment of convective energy transport, I, of 1.0 Hp (pressure scale heights). 


2.1.1. NLTE treatment 


We treat 6706 atomic energy-levels (E-levels) connected by 74550 bound-bound (6 — b) 
transitions of 35 chemical species accounting for various ionization stages of 20 chemical 
elements, including H, He, CNO, and the Fe-group elements that blanket late-type visible 
band stellar spectra, as well as other abundant light metals. We compute the NLTE level 
populations, 'ni(ri 2 ooo )5 and hence the corresponding extinction coefficient, ka('7‘i2ooo), in self- 
consistent multi-level NLTE by solving the system of coupled rate equations of statistical 
equilibrium (SE) consistently with the equation of radiative tr ansfer (RT) in each of the 
relevant bound-bound (5 — b) and bound-free (6 — /) transitions. IShort fc Hauschildtl (120091) 
contains a description of the species treated in NLTE, sources of atomic data, and other 
important details. 


Because NLTE models are more computationally expensive, we only produced NLTE mod¬ 
els and spectra at a subset of our LTE grid, as follows 
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4000 < Teff < 6500 K with AT^s = 500 K, 

-6.0 < [f ] < 0.0 with A[f ] = 1.0, with addition of [f ] = -0.5 


This is sufficient to assess the dependence of NLTE effects on Lick indices throughout 
the grid. Fig. [T]shows NLTE TKinfri^ nno) structures for comparison with those of LTE. NLTE 
radiative equilibrium is co mplex, and lAudersonl fjl989l ) contains a very thorough analysis for 
the case of the Sun, and Short et al. ( 20121) extends the analysis to solar metallicity and 
moderately metal-poor RGB stars. 


NLTE Fe treatment The predicted magnitude of the well-known Fe I NLTE “over¬ 
ionization” , and the resulting predicted brightening of the Fe I-blanketed near-U V a nd blue 
spectr al bands with respect to the rest of the SED (see lShort &: Hauschildtl (120091) and iRutten 
(119861) ) depends on the details of the atomic model of Fe I used in the NLTE Fe treatment. 
More specihcally, the completeness with which high-energy E-levels are included near the 
ionization limit, Xh affects the computed r ate of collisional recombi nation from Fe II, and 
thus the Fe I/Fe II ionization equilibrium (iMashonkina et al. 1120111) . Generally, the more 
E-levels are included for which the atomic energy gap. Ay, between the E-level and xi is 
less than the average collisional energy among particles {kT) throughout the line-forming 
region of the atmosphere, the more accurate the NLTE effect on the computed SED will 
be. For Fe I, xi = 7.9024 eV, and in our 494-level model Fe I atom, the highest lying E- 
level has y = 7.538 eV, for a minimum Ay gap of 0.364 eV. The line-forming region of the 
atmosphere throughout the visible band generally lies at shallower total optical depths, tx, 
than the layer where the continuum value of tx is unity, where T{tx) < Teff. For the warmest 
models in our grid (Teg = 6500 K), kT < 0.560 eV in the line forming region, and we have 
eight E-levels for which Ay < kT, at least in the lower line-forming region. For the coolest 
models in our grid (Teg = 3750 K), kT < 0.323 eV in the line forming region, and we just 
miss having any E-levels for which Ay < kT. We expect our prediction of NLTE Fe I effects 
to be most accurate at the warm end of our grid where collisional recombination into our 
highest E-levels is energetically accessible. At the cool end, the collisional recombination 
rate is artihcially suppressed by the lack of higher-lying E-levels in the model atom. The 
recombination rate is under-estimated for our cooler models, and the NLTE over-ionization 
effect is likely over-estimated. Our NLTE modeled NLTE effects on I values may be thought 
of, cautiously, as upper limits. We plan to expand the PHOENIX NLTE Fe I atom in the 
near future, but this is a significant project in its own right. 























2.2. Synthetic spectra 


Our longer-term goal is to identify useful and overall [^] line diagnostics for high 
spectral resolution from the near UV to the the near IR (NIR). Therefore, we have computed 
synthetic spectra for each of our models for 3000 < A < 26000 A with a spectral sampling, 
AA, set so as to maintain an R value of 300 000 throughout, sufficient to fully resolve 
spectral line cores. This A range includes the NIR J, H, and K photometric bands, in which 
useful line lists of stellar p arameter and abund a nce diagnostics h ave recently been published 
f Bergemann et al. (2012). Cesseti et al. ( 2013 ). Le et al. ( 2011 )1. 


Our synthetic spectra were post-processed by broadening with a Gaussian kernel to R val¬ 
ues of 650 and, following iFranchini et ali (120lOl) . 1800 to match the resolution of the original 
IDS, and SDSS-SEGUE spectroscopy, respectively. A Gaussian is only an approximation to 
the real instrumental spectral profiles of IDS and SDSS-SEGUE spectroscopy, but at this 
stage our study is a differential one to compare the effect of NLTE on Lick indices to that 
of choice of R value as a function of [^]. We do not account for either macro-turbulent or 
rotational broadening. Rotation is expected to be modest in evolved stars of large radius, 
and both effects are expected to be minor at these R values. Fig. [2] shows representative 
synthetic spectra for the models of Fig. [1] convolved to IDS spectral resolution, with the 
Lick indices labeled, and it can be seen that some of the strongest spectral features are still 
significant at [^] values as low as -6.0. Fig. [3] shows the relative flux difference at IDS 
resolution, AF\ = (Fa,nlte — Ea,lte)/Ealte5 models of Tes = 4000 K, and [^] values of 
0.0, -2.0, and -6.0, with the Lick indices labeled. AF\ is generally positive at the A values of 
the Lick indices because most low-y spectral lines from neutral ionization stages of metals 
are weaker in NLTE than in LTE. 


2.3. Synthetic Lick indices 


We use our LTE and NLTE synthetic spectra, convolved to both IDS and SDSS spec¬ 
tral resolution, to compute LTE and NLTE IDS and SDSS Lick indices, /, following the 
prescription of W94. We took the A values defining the latest recommended index and as¬ 
sociated pseudo-continuum bands, similar to the information presented in Table 1 of W94, 
from t he official Lick index W WW site (http://astro.wsu.edu/worthey/html/system.html), 
as did iFranchiui et an (120101) . The IDS indices conform to the we ll-studied Lick index sys¬ 
tem as originally defined, and are directly comparable to those of iGorgas et al\ (119931) and 
W94, and serve as a check on our procedures, as well as allowing us to asse ss the impact of 


NLTE effects at IDS resolution. The SDSS indices are comparable to those of iFranchini et al. 
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fl2010l) and allow an assessment of NLTE effects at somewhat higher resolution typical of 
more modern spectroscopic surveys such as SDSS and LAMOST. Furthermore, comparing 
the LTE IDS and SDSS indices allows an assessment of the dependence of the I sensitivity 


to Tefj and [^] on R value. 


Following iGorgas et all fjl99311 and W94, we also computed the photometric V — K index 
for the models as an observational surrogate for the independent parameter Teg. V — K 


has been found to be relatively insensitive 

to log^f and [^], and a good proxy for over 

). For consistency with W94, we use the V- and 

the GK star range ( 

Bell & Gustafsson 

1989 

A-band hlter dehnitions of 

Johnson et aL ( 

1966 

) and calibrate the index with a single-point 


calibration of the 


[ii 


= 0 . 0 mod els at T^ff = 4000 K to the ( V — K) — Teg relation given in 


Table 4 of iRidgwav et all fll980l) . The iRidgwav et al\ (1198011 {V — K) — Teg relation is for 
giant stars, and a Teg value of 4000 K is near the center of their calibrated range, and 
overlaps with the Tes range of our grid. We always use the LTE model V — K color on the 
grounds that it is serving as an independent variable in this analysis, and the LTE grid is 
more complete. 


It is worth reiterating remarks made by previous investigators about the particular diag¬ 
nostic utility of those indices that are expected to be most useful in this investigation: 


Fe4383 and Fe4668 (W94) and Fe5270 and Fe5335 ( Gorgas et ai?.lll993l ) are [^]-sensitive 
with a range in I value significantly greater than measurement uncertainty, and are expected 
to be especially useful here if they remain detectable down to XMP metallicities (note that 
Fe4668 has signihcant contributions from Mg, Gr, Ti, and G 2 ). Ga4227 really is dominated 
by Ga (whereas Ga4455 is more influenced by Fe-group lines), thus providing one of the few 
atomic indices not heavily affected by Fe, and is somewhat sensitive to overall [^] (W94) as 


well as [§1 given that Ga is an a-process element. GN 2 is a modification of GNi designed 


LHJ 


to avoid contamination from the Hh line, and is strongly dependent on overall [^] for giants 
(but less so for dwarfs) (W 94). By contrast to the preceding, II/3 has been found to be most 
strongly dependent on Teg (IGorgas et a/.lll993fl . and thus provides a val uable complernentary 
diagnostic. Indices Mgi (dominated by MgH) and Mg b were found by IGorgas et al\ (jl993l ) 
to be usefully sensitive to log^f, and we include them in our investigation to re-assess their 
jj j sensitivity (note that Mg 2 includes contributions from both MgH and the Mg b 


Rff and 


lines, so is less pure a signal of either). Na D is known to be signi hcantly contaminat ed by 
interstellar (ISM) extinction, which complicates its interpretation (IGorgas et a/.lll993l) . 
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Franchini et al\ 020101) investigated the influence of enhanced a-element abundances on 


model I values. They found that for Tgfj > 4250 K the most a-sensitive indices are CNi, 
CN 2 , CaHK, Ca4227, Fe4668, Mgi, Mg 2 , and Mg b, and that the least a-sensitive indices 
are G4300, Ca4455, Fe4531, Fe5015, Fe5782, and H/1. For < 4250 K the situation seems 
more complex, but CaHK, Ca4227 and Mg b remain among the most a-sensitive, and G4300 
and Ca4455 remain among the least a-sensitive indices. 


3. Results 

We caution that because our models have scaled-solar abundances, in the discussion 
that follows the [^] parameter is effectively identical to the [^] parameter in the internal 
context of our modeling and analysis, whereas in a-enhanced metal-poor RGB stars with 
non-solar abundance distributions, [^] differs from [^]. This distinction is expected to be 
most important for those Lick indices that are Mg- or Ca- dominated, and less so for those 
that are Fe-dominated. 


3.1. K - K{Tes) and K - R:([f ]) relations 


Fig. m shows the LTE and NLTE model V — iF(Teff) relation for our range o f mode l 


[^] values, over-plotted with the V — K{Tes) relation for giants of iRidgwav et al\ (Il980l) . 


The model V — K{Tes) relation flattens with decreasing [jj] value, which is to be expected 


rMi 


because line blanketing extinction in the V band increases more rapidly with increasing jj j 
value than does that in the K band. Within the range of overla p in Tpff (3750 t o 5000 K), our 
model V — K(Tes) relation at [ ^] = 0.0 close l y trac ks that of iRidgwav et al\ fll980l) . but is 


slightly steeper. However, the iRidgwav et al\ fjl980l) sample of red giants probably includes 


stars of [^] < 0.0, so it should have a flatter V — iF(Tefj) relation than that of [^] = 0.0. 
NLTE effects are negligible, which is to be expected given that these are broad-band colors 
that average the effects of many spectral lines, and that line blanketing opacity is already 
considerably reduced in the V band as compared to the B and U bands. 


3.2. XMP indices 

Table[2]displays, for a selection of Tgfj values spanning our grid, the range of [^] values for 
which each index, /, is a sensitive [^] indicator as judged by the criterion that A[^] x > 

(Tworthey, where A[^] 1 and Uworthey is an observational uncertainty described below. Lick 
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indices that meet this criterion and remain strong enough at VMP-to-XMP metallicities 
([^] < —4.0) to be detectable are considered to be “Lick-XMP” indices. Generally, all the 
“metallic” atomic and molecular indices have values that increase with decreasing T^s 

value, and we expect that Lick-XMP indices will be more readily identihable at the cool end 
of our grid. 


Table [2] includes indications for indices that become “pathological” over a signihcant [^] 
range at any values by becoming multi-valued, or negative in the case of those indices 
that are in linear W\ units. These pathologies often reflect complications in the bracketing 
pseudo-continuum wave-bands as dehned by the Lick index standard rather than with the 
central feature itself, and most often appear at higher [^] values. 


At Teff = 3750 K we have found hve Lick indices that meet our Lick-XMP criterion all 
the way down to [^] = —6.0 (Fe5270, Fe5335, Fe5406, Mg b, Na D), and four others that 
meet the criterion down to [^] = —5.0 (Fe4383, Fe4531, Fe5015, Mgi). At Tes = 4500 K 
there are none that meet the criterion down to [^] = —6.0, and four of the above nine that 
still meet the criterion down to [^] = —5.0 (Fe5335, Fe5406, Mg b, Na D). These nine 
Lick-XMP indices are indicated in Figs. [2] and [3] with identihcation labels that are off-set 
from the rest. By Tes = 5000 K (and hotter) there are no indices that meet our Lick-XMP 
criterion. (Furthermore, we caution that Mg b becomes multi-valued as a function of [^] for 


T'efr > 5000 K.) The presen ce of Fe4383, Fe5270 and Fe5335 among our Lick-XMP indices 
is not surprising given that iGorgas et al\ f 19931) identihed them as strong [^]-indicators. 
Na H is an interesting member of our Lick-XMP indices in that it is not Fe-dominated, 
but we caution, again, that its usefulness is compromised by signihcant ISM extinction. As 
noted in Section [21 our treatment of Mg b and Mgi is least accurate because we neglect 
a-enhancement in the current investigation, and Mg is an a-process element. Therefore, the 
following results for these two indices are most suspect and require further investigation. 


Figs. |5]to[8]show the modeled /([^]) relation at = 4000 K, and Figs. 191 to IT2] show the 
modeled I{Teg) relation at [^] = 0.0 for four of our Lick-XMP indices, Fe5270, Fe4383, Mg b, 
and Na D. The latter are the only two Lick-XMP indices not dominated by Fe. Results are 
shown for spectra computed with the R values of SDSS and IDS spectroscopy, as computed 
in LTE and NLTE. For comparison, we have over-plotted observati onally derived I(Tpfr) 
and_/([^]) points for a wide range of cool giants from the catalog of lWorthev fc Ottaviani 


fll997l) with catalog logg values between 1.0 and 3.0, and and 
±1.0, respectively, of the plotted models. 


HI 


within ±1 0 0 K a nd 
Note that the IWorthey fc Ottaviani (1997) I 
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da ta only includes stars of [^ ] > —1.0, and that no calibration of onr / valnes to those 
of IWorthey fc Ottavianil (jl997|) has been performed. Nevertheless, from Figs. [7] and El 
the agreement between onr modeled Mg b and Na. D I values and the measurements of 


Worthey fc Ottavianil (119971 ) for stars of Tes = 4000 ± 100 K is assuring. Figs. [5] to El show 


that these indices generally satisfy our Lick-XMP criteria for ^eff values in the cool part of 
our grid. Fig. [TOl shows that the utility of Fe4383 is somewhat compromised - it is double¬ 
valued as a function of Tes at the cool edge of our grid. However, it is a sensitive and useful 


[^] discriminator otherwise, and we include it because there are few indices that qualify as 
Lick-XMP at all. Figs. 191 to fT^ show that, as discussed above, Lick-XMP indices quickly 
become weaker and lose their ability to discriminate among [^] values for ^ 5000 K. Figs. 
[5] through [12] also include an indication of the “observational uncertainty” as determined by 
W94, cTworthey (see below), to aid in assessing the signihcance of AJ differences. 


The effect of NLTE is complex in that every Lick index is really a compound feature 
caused by signihcant spectral lines from several species. Moreover, although the bracketing 
pseudo-continua used to dehne the indices were chosen to be relatively insensitive to stellar 
parameter values, NLTE effects on line strengths in the bracketing regions will, in principle, 
also play a role in the overall effect of NLTE on the computed I value. However, Fe5270 is 
typical of our results for the Fe-dominated indices in tha t the effect of the well-known NLTE 
over-ionization of Fe I in late-type stars flRuttenl Il986l) leads to smaller I values at every 
Teff — [^] combination. Figs. [13] and [H] show the size of the NLTE effect, AJ = /nlte — 7lte, 
as a function of for Fe5270 and Mg b. The effect of NTLE on strong low x Fe I lines is to 
weaken them (a negative correction to modeled I value) as a result of NLTE over-ionization. 
In late-type stars the effect will generally be maximal where the discrepancy, AT, throughout 
the line forming region between the radiation temperature, Tr sh ? of the photo-ionizing near- 
UV band radiation in the NLTE treatment, and the local kinetic temperature, Txin, that 
determines the ionization balance in the LTE treatment, is largest. For Fe I in giants, the 
discrepancy is largest around 5000 K, and decreases in magnitude for both lower and 


higher Teg value (see iRutten] (jl986[) for a thorough analysis for the case of Fe in the Sun). 


For Mg I (Index Mg b) the NLTE correction to the modeled I value is also negative, 
but increase i n magnitude w i th de creasing Tgfj throughout this Teg range, especially for 
Teff < 4000 K. lOsorio et al. fl2015l) very recently conducted a thorough NLTE analysis of 
Mg I in late-type stellar atmospheres, including an investigation of H I collisional cross- 
sections and electron-exchange reactions. For the A 5184 line, they found that NLTE effects 
lead to a weakening of the modeled line (hence the positive abundance correction in their 
Fig. 10), by an amount that depends on choice of atomic data, but can be as much as 0.4 
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dex at Teff = 4500 K, log^f = 1.0, and [^] = 0.0, which is qualitatively consistent with our 
result. 


Interestingly, the magnitude of NLTE effect on the computed / value is comparable to 
that caused by changing spectral resolution {R value of IDS vs SDSS). Na D is an exception 
because it is so broad that it is minimally affected by the choice of R. Computing I from 
higher i?-value spectral (he. that of SDSS) can either increase or decrease the / value, 
depending on the index. The same remarks as made when considering NLTE effects above 
also apply: namely, that the effect of R value on any given index will depend on how the 
bracketing pseudo-continua are affected as well as the central feature itself. Altogether, the 
vertical spread in I values at each abscissa can be taken as an approximate indication of 
“spectroscopic and modeling physics uncertainty”. 


3.3. Polynomial fits 

Based on the approach taken by W94 with observed IDS spectra and observationally 
determined stellar parameters, we have used multi-variate linear regression to determine 
ten polynomial htting (regression) coefficients (or model parameters), {Cn}, for each index, 
/, for a polynomial htting function, P 3 , that accounts for all terms, including cross-products, 
up to third order in [^] and log^, where 9 = 5040/Teff, 


I PS P;^ = Cq + C'liTy] + ^2 log 9 + Cs)—-l- C 4 log^ 9 + C'slTy] ^ 

rl rl rl 

+ C 7 log^ 9 + log 9 + ^ 9 [^] log^ 9 (1) 

Note that, as per convention, the units of I for the molecular indices CNi, CN 2 , Mgi, 
Mg 2 , TiOi and Ti 02 are magnitudes, and those for the remaining indices are A. Moreover, 
following W94, in the special case of the TiOi and Ti 02 indices, which exhibit a rapid 
increase in strength with decreasing Teg near the lower limit of our Teg range, we ht Eq. [T]to 
log / so that I is being ht by exp P^. The htting is performed with the intrinsic REGRESS 
procedure in Interactive Data Language (IDL) and is achieved by minimizing the hgure 
of merit assuming that the uncertainty, cri(/), of all “data” points R, is unity. We note that 
IDL installations include the source code for all intrinsic procedures, and we have been able 
to critically inspect the REGRESS procedure. We use REGRESS to compute nine linear 
regression coefficients and a constant term in Equation [T] for nine basis functions that consist 


14 


of the powers of the independent parameters (Teg, [^]) and their products up to third order. 
This amounts to htting a model with nine parameters to 90 data points (90 computed I values 
for ten Tgg and nine [^] values), thus having 81 degrees of freedom. This is consistent with 
the htting method of W94 - we note that their Equation 4 appears to be the standard formula 
for because their summed square deviations are weighted by their inverse observational 
uncertainty, 1/cr^, although they have labeled their hgure of merit “rms^”. We note that 
because we have no proper data uncertainties (he. “measurement errors” in data modeling), 
(Ji(/), we cannot properly propagate errors to compute uncertainties for the htted model 
parameters, Cn - rather we compute “htting uncertainties”, a, post hoc from the calculated 
value of cind the number of degrees of freedom. Assuming these “htting uncertainties” 
rehect errors that are normally distributed, they may be interpreted as 68% conhdence 
intervals for the htted values of the corresponding C^- 


Tables [3] and m present these {Cn} values for LTE spectra of IDS and SDSS resolution, 
respectively, for the nine indices that were identihed above as good Lick-XMP diagnostics, 
in the same format and numerical precision as that of Table 2 of W94 (“Data for stars of 
3570 < Teff < 5160 K”) for direct comparison. Table [5] presents the results for IDS spectral 
resolution of performing the same third order multiple linear regression with the model 
V — K color in place of log 6*. We also present the values of the standard deviations, cr, 
computed for each Cn parameter from and the value of the reduced given 81 degrees 
of freedom, although we caution that in the absence of proper measurement errors, y^ is 
not really a goodness-of-ht hgure of merit. Figs. [15] through [18] show the comparison of the 
polynomial hts to the modeled /(Teg) relation for [^] = 0.0 at SDSS and IDS resolution, and 
the residuals. For SDSS resolution, we also show both the htted relation and the residuals 
computed with the 1 -a “htting uncertainties” (see above) added and subtracted from each 
of the Cn values to illustrate these two limiting cases. Figs. [19] through [22] present the same 
information for the polynomial hts to the modeled /([^]) relation for = 4000 K. 


There are a number of important diherences between our approach and that of W94: 


a) Our On values complement those of W94 by being based on spectra of SDSS resolu¬ 
tion rather than those of IDS resol ution, and are ni o re rel evant to both SDSS and LAMOST 
spectra and to the investigation of mranchini et al\ fj2010l) . 


b) log g is not a htting parameter, so we do not have cross-product terms that capture 
the dependence of I on the product of log^f or any of its powers and other parameters or 
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their powers. 


c) Table 2 of W94 contains coefficients for fits in the Tgfj range of 3570 to 5160 K, 
whereas onr fits apply to the range 3750 to 6500 K. The difference at the high end 
is necessary for us to accommodate the Teg range of interest for detected halo red giants. 
Because the lower limit of our Tgg range is signihcantly higher than that of W94, results for 
our TiO indices are especially suspect, and not comparable to W94. W94 included stars of 
Tefj > 5160 K in their hts for warm and hot stars (5040 to 13260 K, their Table 3). This is well 
beyond the limits of our red giant grid and we are not able to compare to their “hot star” hts. 

d) W94 carry out a careful statistical F-value test of the goodness-of-ht to separately 
determine whether the addition of each successive Cn term in Eq. [T] led to a statistically 
signihcant change in the htted I value for each index. As a result, many of the values 
in their Table 2 are blank because, presumably, including the corresponding term in the 
htted function led to an insignihcant reduction in the variance. We have chosen to simply 
let the regression hnd values for all ten terms in Eq. [T] consistently for all indices, with 
the expectation that terms that are of low signihcance will have small htted values of the 
corresponding C^- Our expectation is that a third order ht is of low enough order that we do 
not expect high order spurious solutions to compromise the ht, and our situation is simpler 
than that of W94 in that logg is not a htting parameter. 


Error analysis W94 describes the uncertainty in determining an I value as a “typical 
rms error per observation”, which we denote aworthey For reference, we have included an in¬ 
dicator of the magnitude of cJworthey for each index in Figs. [51 through fT^ W94 quantihes the 
uncertainties in the htted I values with a “residual” rms value in units of the observational 
uncertainty, cxworthey We are not working with observational data, and we quantify the un¬ 
certainties in our hts, a, with the quadrature sum of the 1 -a uncertainty estimates computed 
for each value from the multiple linear regression procedure (as described above), and 
have included them in Tables [3] through [5l Generally, the hrst order coefficients, Ci for [^] 
and (7*2 for Teg, from the ht to SDSS resolution spectra are larger than those form the ht 
to the IDS resolution spectra. This is to be expected because in spectra of higher R value 
the hrst order dependence of the strength of spectral features on stellar parameters is less 
diluted by the re-shuffiing of information within the instrumental spectral prohle. 


The coefficient of the log^ 9 term, Gy, generally had, by far, the largest 1-a uncertainty 
value of all the Gn values, and dominates our quadrature sum a values. Generally, the 
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parameter with next largest 1-a uncertainty was the log^ 9 term, C 4 , but it was much smaller 
than that of Cj. Table [3] shows that the magnitude of Cj is generally larger than any of the 
other coefficients, and this is consistent with what was reported in Table 2 of W94. We also 
found that the value of Cj was the most sensitive to spectral resolution, differing by as much 
as a factor of four between the fits to IDS and SDSS resolution spectra. We conclude that 
the terms in non-linear powers of log^ were generally least well fit. However, for these GK 
stars, 6 = 5040/Tefj is of order unity, and the squares and cubes of the independent variable 
log 6 * is much less than unity. We have found that the 1-a uncertainty of the C 4 and Cj 
terms for the analogous V — and V — terms from the hts with V — K m lieu of Te^ 
are much smaller and are consistent with those of the other coefficients, and from Table 
|5] it can be seen that the corresponding total a values are smaller. 


For our nine XMP indices, we compare our htted Gn values for the zeroth and hrst order 
terms of the £t to IDS resolution spectra with those of W94, given the four caveats listed 
above. For the hts with log 6 ' as an independent parameter, comparable to Table 2 of W94, 
for hve of six indices designated “Fe”, our htted Cq value is consistently smaller than that 

of W94, and ranges from about 0.6 to 0.75 of the W94 value, and for Fe5270 we are in 

very close agreement. For Na D, our Cq value is also smaller, but is in closer agreement 
with that of W94. Mgi is our only index for which our Cq is larger than W94, by about a 
factor of 1.5. For Mg b we hnd a negative value of Cq where W94 finds a positive value. In 
both cases, the magnitude of Cq is about unity. For the log^ coefficient, C 2 , our values for 

Fe4531, Fe5335, and Fe5406 are close to the values of W94. The remaining C 2 values are 

also generally within a factor of two, greater or less than, and of the same sign as those of 
W94. The two exceptions are Na D for which our C 2 value is about three times larger, and 
Mgi for which our C 2 value is positive and that of W94 is negative, although the magnitudes 
are with a factor of two. For the [^] coefficient, Ci, our values for Fe5335, Fe5406, Mg b, 
and Mgi agree closely with those of W94. For the remaining XMP-Lick indices, including 
Na D, our Ci values diher by as much as a factor of 2.5, greater or less then, those of W94. 


For the hts with V — K as an independent parameter (in lieu of log 6 '), only four of our 
nine Lick-XMP indices appear in the comparable Table 4 of W94 (Fe4383, Fe4531, Fe5015, 
and Fe5406). For Fe5015, W94 does not present a Cq value, but our C 2 ( 7 (H — K)) agrees 
closely with theirs. For Fe4531, W94 does not present a C 2 value, and our Cq value is about a 
factor of two larger. For both of these indices, our Ci values (7([^]), at constant V — K this 
time) are close to those of W94. For Fe4383 and Fe5406 the situation is disconcerting and 
puzzling - we hnd both Cq and C 2 values with the opposite sign and a diherence in magnitude 
ranging from two to four. For both of these indices, W94 have no Ci value for the term in 
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[^], indicating, presumably, they found insignificant /mear dependence of / on [^] at fixed 
V — K, and, consistently, we find modest Ci values of -0.4489 and -0.0986, respectively. 
We have been unable to identify the reason for the gross discrepancy in htted polynomial 
coefficients for Fe4383 and Fe5406, beyond those expressed in the caveats itemized above. 


3.4. Special indices 


CaHK ISevern. Worthev &: Brilev! (20051) introduced a new Lick-type index, which they 


designate CaHK, and iFranchini et al\ (120101) modihed the dehnition by changing the blue 
pseudo-continuum band-pass and removing the central 5 A around the line cores to remove 
potential chromospheric emission. IFranchini et al\ (120101) found that CaHK had the advan¬ 
tage of being sensitive to a-enhancement, as well as being less influenced by Fe than most of 
the other atomic indices. The index seems to shows promise as a Lick-XMP [^] diagnostic 
because the line is very strong. We do not remove the central 5 A from our index compu¬ 
tation because our models have outer atmospheres that are in radiative equilibrium and do 
not have chromospheric emission. We have conhrme d that we can reproduc e the qualitative 
double-valued behavior of I{6) exhibited in Fig. 7 of iFranchini et al\ (120101) . We have found 
that the index is a good Lick-XMP diagnostic for our models at the cool edge of our grid, of 
Teff value equal to 3750 to 4000 K. However, unfortunately, CaHK is either double-valued as 
a function of [^] for Teff in the range of 4250 to 5500 K, or becomes negligible and insensitive 
to [^] for [^] < —4 in the Teff range of 5500 to 6500 K. We can only recommend CaHK as 
a Lick-XMP index for Teff < 4000 K. Because the behavior of CaHK is double valued as a 
function of both Teff and [^] throughout much of the ^eff — [ plane, the and partial 
derivative values are probably not useful for interpolation, and we have not included them. 


TiO We have included the TiOi and Ti02 indices in our analysis because they are part 
of the Lick system, but they are only signihcant in strength and in [^]-sensitivity for the 
coolest part of our grid, and then only at the highest [^] values. Ti 02 is the stronger of the 
two indices and remains sensitive to [^] in the -2.5 to 0.0 range to Teff values as high as 4250 
K. Neither TiO index has pathologies such as being double-valued. Note that we £t Eq. [T] 
to log / to account for the strong Tgff and [^] dependence of TiOi and Ti02, following W94. 
Because our grid does not extend to the low Teff values included in the W94 fit, where TiO 
is strong, our and partial derivative values are not comparable to W94, and should be 
treated with more caution that those of the other indices, and we do not include them here. 
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3.5. Partial derivatives with respect to V — K and [^] 

We have used our {Cn} values to compute the partial derivatives 100 K x g|^|[M/H], 
^ ^ a(v-K) l[WH]) and 0.5 x ap^l(v-ir) for all Lick indices modeled 

in LTE, in both index units and in units of a, as dehned above, for IDS and SDSS resolution 
spectra. For example, the partial derivative with respect to log 6* at constant [M/H] can be 
found from 


g|;D|,„/„,=C2 + 2C4log» + C,|M| + 

3C7 log^ 9 + + 2C9 log f?[—]. 

Then, the partial derivative with respect to T^e follows from 


( 2 ) 


dl 


dl 

dlogO 


log e/Teff. 


(3) 


Tables |6] and [7] present the values for our nine identified Lick-XMP indices for SDSS 
resolution only, and are comparable to Tables 7A and 7B of W94. The values in Table |6l 
provide an indication of by how much the measured value of I will differ between two stars 
that differ by ATefj ~ 100 K at each value of [^], and by how much / will differ between 
two stars that differ by A[^] 0.5 at each value of Teg. Alternately, these derivatives can 

be used to estimate the change in inferred T^e or [^] as a result of the change in computed / 
caused by NLTE effects (he. where AJ = Jnlte — 7lte)- The partial derivative values that 
are in units of a give an indication of the significance, or detectability, of these changes in I 
value given the formal uncertainty AJ a. 


As an example of the implications of NLTE effects for [^] determination, from Table |6l 
for Fe5270 the quantity 0.5 x l^ff ~ 1-2717 at (Tgff, [M/H]) = (5000 K,0.0), and from 
Fig. [IHl the computed change in I caused by NLTE effects, AJ, at 5000 K is ~ —0.05. This 
corresponds to an LTE model [M/H] value that is smaller by 

A[M/H] ^ 0.5 X -0.05/1.2717 ^ -0.02 (4) 


Inversely, fitting a given observed / value with NLTE as compared to LTE 
models would require a compensating model value of [M/H] that is ~ 0.02 larger, 
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consistent with the sign of the change in inferred [^] at fixed / value for = 4000 
K shown in Fig. [5l We emphasize that this is an estimate of A[M/H] based on the 
modeled LTE value of ^[m/h] l^eg • ^ more accnrate estimate wonld follow from a NLTE valne 
of g[M/H] importance of this consideration depends on the magnitnde of the 

difference between ^[m/h] Ihsg valnes compnted from NLTE and LTE model grids. However, 
compntation of NLTE partial derivatives that are comparable to those of LTE reqnires a 
NLTE model grid that inclndes all the same (Tgff, [M/H]) points as the LTE grid. Cnrrently, 
onr NLTE model grid covers only a snbset becanse of the larger compntational cost of NLTE 
modeling, and its only pnrpose is to spot check the effect of NLTE on compnted I valnes at 
select grid point. 


Conclusions 


Using a grid of red giant synthetic spectra that extends from solar- to XMP-metallicity 
([w] “ —6.0) we have identified nine of the original 21 Lick indices, designated Lick-XMP 
indices, that remain significantly detectable and significantly sensitive to [^1 down to XMP 


■ hj 


valnes (at least [^] = —5.0) for giants of T^s < 4500 K. For warmer late-type giants, all Lick 
indices become nndetectable or insignificantly sensitive to [^] before [^] decreases to -4.0. 
The Lick-XMP indices shonld be the most nsefnl ones for characterizing very old “fossil” 
stars that formed very early in the history of the Galax y, and in other galaxies. We also 
investigated a new ^er Lick-type ind e x, Ca HK, introdnced bv ISevern, Worthev fc BrilevI (1200511 
and developed by iFranchini et ali (120101) as a potential Lick-XMP index, given its strength. 


However, for CaHK, / is donble valned as a fnnction of [^] and over mnch of onr grid 
and its nsefnlness is restricted to the cool edge of onr grid (Tgff < 4000 K). 


For onr LTE grid of SDSS resolntion spectra we present polynomial coefficients, {Gn}, to 
third order in the independent variable pairs (log6', [^]) and ((U — K), [^]) derived from 
mnlti-variate linear regression, approximately comparable to those of W94 for IDS reso¬ 
lntion spectra. We present the partial derivatives ^[m/h] WeS, I[m/h]; and 

a[M/H] I (v-K) compnted from onr {Gn} valnes. 

For Fe-dominated Lick indices, the effect of NLTE is to generally weaken the valne of I 
at any give [^]) valne. To pnt the magnitnde of the NLTE effect into context, the 
change, AJ, cansed by NLTE effects is generally comparable to the change that resnlts from 
compnting I from spectra of SDSS resolntion rather than that of IDS. The partial derivatives 
can be nsed to estimate a change in inferred T^g at fixed [^], or a change in inferred [^] at 
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fixed Teff, resulting from a change in / {e.g. as caused by NLTE effects) at any (Tgff, [^]) 
value pair throughout our grid. For example, from Figs. [9] and [10], for stars of inferred 
Teff > 4200 K (6* < 1.2), an Fe-dominated / value computed in LTF that is too strong might 
be compensated for by inferring a Teg value that is too large, for hxed inferred [^]. 


CIS is grateful for NSFRC Discovery Program grant RGPIN-2014-03979. The calcula¬ 
tions were performed with the facilities of the Atlantic Computational Excellence Network 
(ACEnet). 
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Fig. 1.— TKin(Ti 2 ooo) structures for models of T^s = 4000 K and [^] values of -0.5, -2.0, 
-4.0, and -6.0. Gray-scale: darker line indicates larger [^] value. Results are shown for LTE 
(dashed line) and NLTE (solid line) models. 
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Fig. 2.— Synthetic spectra for the LTE models of Fig. [H convolved to IDS spectral resolution 
for clarity. The gray-scale is as in Fig. [1] Index labels that are offset upward refer to our 
nine Lick-XMP indices. 
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Fig. 3.— Relative difference of synthetic spectra convolved to IDS resolution showing the 
comparison of NLTE to LTE results for models of Teg = 4000 K and [^] values of -0.5, -2.0 
and -6.0. The gray-scale and index label positions are as in Fig. [2l 
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I 



4000 5000 6000 7000 


Fig. 4.— i/ —iF(Teff) relation. Synthetic colors from LTE (solid gray-scale li nes and crosses) 


and N LTE (dashed gray-scale lines and diamonds), and the calibration of 
f lQSnh (solid black line and triangles). The synthetic colors were tied to the 
f 198ol) relation with a single-point calibration of the [w] = 0.0 models at = 4000 K 


Ridgwav et al. 


Ridgway et al 


(vertical line). Gray-scale: Darker line indicates larger [^] value throughout model grid 
range of -6.0 to 0.0. Note that we only compute NLTE I values at a subset of the LTE grid 
- see text. 
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Fe5270 ^,^^=4000 K 



Fig. 5.— Index Fe5270, one of our identified Lick-XMP indices: / as a function of [^] 
at Teff = 4000 K. Vertical line in middle left: Indication of the observational uncertainty, 
cworthey, of W94 - See text. / has been computed from spectra with R values of SDSS (solid 
line) and IDS (dashed line) spectroscopy that were computed in LTE (triangles, squares) 
and NLTE (crosses, diamonds) - see text. The vertical spread in / values at a given [^] value 
can be taken as an estimate of “spectroscopic and modeling physics uncertainty”. Note that 
we only compute NLTE I values at a subset of the LTE grid - see text. For com parison we 
show the observationally derived / and Teg values of lWorthev &: Ottavianil (1199711 for giants 
(1.0 < l ogq < 3.0) of Tpff with i n ±1 00 K (“X” symbols). No calibration of our I values to 
those of lWorthev fc Ottavianil (119971) has been performed. 
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Fe4583 ^,^^=4000 K 



Fig. 6.— Same as Fig. |5l but for index Fe4383. 
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Mg_b r,ff=4000 K 



Fig. 7.— Same as Fig. |5l but for index Mg b, an example of one of our identified Lick-XMP 
indices, not dominated by Fe. 
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Na_D r,ff=4000 K 



Fig. 8.— Same as Fig. |5l but for index Na D, another example of one of our identified 
Lick-XMP indices, not dominated by Fe. 
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0.8 1.0 1.2 1.4 

d = 5040/rEff 


Fig. 9.— Same as Fig. [5l but for J as a function oi 9 = 5040/ Tpff at [^] = 0.0. For 
compa rison we show the observationally derived / and [^] values of lWorthev fc Ottaviani 
Jl997h for stars of [f] within ±1.0 (“X” symbols). 








/ (A) 
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0.8 1.0 1.2 1.4 

0 = 5040/rEff 


Fig. 10.— Same as Fig. |9l but for index Fe4383. / is double valued at the cool edge of the 
grid, making its use there problematic, but is generally useful as an [^] diagnostic otherwise. 
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0.8 1.0 1.2 1.4 

6 = 5040/rEff 


Fig. 11.— Same as Fig. [9l but for index Mg b. 
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Na_D [M/H] = -0.0 



Fig. 12.— Same as Fig. [9l but for index Na D. 




A Index (A) 



Fig. 13.— Index Fe5270: The difference between / computed in NLTE and that in LTE, 
AJ = /nlte —-^LTE, as a function of 6* = for [^] = 0.0. Results for / derived from spectra 
of SDSS (crosses) and IDS (diamonds) resolution, R. The horizontal dotted line indicates a 
AJ value of zero. 




A Index (A) 
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Mg_b [M/H] = -0.0 



Fig. 14.— Same as Fig. [131 except for index Mg b. 
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Fe5270 [M/H] = -0.0 



Fig. 15.— Fe5270: I{9) based on multiple linear regression fits for /(log 6*) at [^] = 0.0. 
Upper panel; Fitted relation (solid line) to modeled / values (triangles) at SDSS resolution, 
and similarly at IDS resolution (dashed line and squares). The dotted lines show the htted 
relation for SDSS resolution computed with htting coefficients, C^, that have had their 1-a 
error from the htting procedure (see text) added and subtracted from them. Lower panel: 
The residual values for SDSS and IDS resolution (solid and dashed lines, respectively). The 
dotted lines are the residuals for SDSS resolution for a htted relation computed with 
values with their l-a errors added and subtracted (see upper panel caption). The horizontal 
gray line indicates a residual of zero. 
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Fe4383 [M/H] = -0.0 



1.0 

< 

X 0.5 

0.0 

I 

-0.5 
- 1.0 
0 . 



Fig. 16.— Same as Fig. [151 but for Fe4383. 
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Mg_b [M/H] = -0.0 




e = 5040/7’.,, 


Fig. 17.— Same as Fig. [151 but for Mg b. 
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No_D [M/H] = -0.0 

T-1-r 




0.72 0.84 1.01 1.26 

B = 5040 / 7 ’.,, 


Same as Fig. [151 but for Na D. 
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Fe5270 7’e„=4000 K 




-6 -5 -4 -3 -2 -1 0 

[M/H] 


Fig. 19.— Same as Fig. [151 but for the /([^]) relation at = 4000 K. 









Mod 


41 


Fe4383 7’e„=4000 K 




Fig. 20.— Same as Fig. [191 but for Fe4383. 
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Mg_b 7’e„=4000 K 




Fig. 21.— Same as Fig. [191 but for Mg b. 







Mod 
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-6 -5 -4 -3 -2 -1 0 



Fig. 22.— Same as Fig. [191 but for Na D. 
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Table 1. Index values, J, at SDSS spectral resolution from LTE models for the nine 
Lick-XMP indices (see text) for models of = 4000 K and logg = 2.0, in which the I 
values are especially large, and select [^] values computed with scaled solar abundances 
(Erst number in each data pair) and with a relative enhancement of a-element abundances 
by +0.4 (second number in each data pair). We note that for the scaled-solar abundance 
models, [^] is identical to our [^] grid parameter. 


Lick index, / 


[f] 

Fe4383 

Fe4531 

Fe5015 

Mgi 

Mg b 

-2.0 

6.235, 5.454 

4.041, 4.302 

5.167, 5.257 

0.189, 0.275 

1.845, 2.449 

-4.0 

1.098, 1.002 

0.520, 0.716 

0.647, 0.891 

0.006, 0.016 

0.433, 0.670 

[f] 

Fe5270 

Fe5335 

Fe5406 

Na D 


-2.0 

3.040, 2.988 

2.635, 2.523 

1.923, 1.797 

1.217, 1.054 


-4.0 

0.488, 0.545 

0.500, 0.546 

0.374, 0.404 

0.271, 0.288 
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Table 2. Range of [^] value for which Lick index, /, is a significant [^] diagnostic for 
differences of A[^] 1 for select Tes values, and corresponding {V — K) values at 

[A] = 0 . 0 , at SDSS spectral resolution. 


Index 

3750 

4.33 

TeS 

4500 

V -K at 

2.58 

(K) 

5000 

[M] = 0.0 
2.02 

6500 

0.89 

CN2 

1 

1 

o 

-1.5, 0.0 

-0.5, 0.0 

is 

Ca4227 

-4.0, 0.0 

mv, neg 

mv, neg 

-3.0, 0.0 

G4300 

mv 

-3.0, -5.0 

-1.5, -4.0 

0.0, -3.0 

Fe4383* 

-5.0, -1.0 

-4.0, 0.0 

1 

1 

o 

-3.0, 0.0 

Ga4455 

-4.0, 0.0 

-3.0, 0.0 

-3.0, 0.0 

-2.0, 0.0 

Fe4531* 

-5.0, 0.0 

-4.0, 0.0 

-3.0, 0.0 

-3.0, 0.0 

Fe4668 

-1.0, 0.0 

-1.5, 0.0 

-1.5, 0.0 

-1.0, 0.0 

Fe5015* 

-5.0, 0.0 

-4.0, 0.0 

-4.0, 0.0 

-3.0, 0.0 

MgR 

1 

1 

o 

-3.0, 0.0 

-2.0, 0.0 

-1.0, 0.0 

Mg b* 

-6.0, 0.0 

-5.0, 0.0 

mv 

mv 

Fe5270* 

-6.0, 0.0 

-4.0, 0.0 

-3.0, 0.0 

-2.0, 0.0 

Fe5335* 

-6.0, 0.0 

-5.0, 0.0 

-4.0, 0.0 

-3.0, 0.0 

Fe5406* 

-6.0, 0.0 

-5.0, 0.0 

-3.0, 0.0 

-2.0, 0.0 

Fe5709 

-4.0, 0.0 

-3.0, 0.0 

-3.0, 0.0 

-1.0, 0.0 

Fe5782 

-3.0, -1.0 

-2.0, 0.0 

-2.0, 0.0 

-1.0, 0.0 

Na D* 

-6.0, 0.0 

-5.0, 0.0 

-4.0, 0.0 

is 

TiOi 

-3.0, 0.0 

-1.0, 0.0 

-1.0, 0.0 

is 

TiOs 

-4.0, 0.0 

-3.0, 0.0 

-2.0, 0.0 

-1.0, 0.0 


is, mv, neg: I suffers from one or more of several pathologies over a signihcant range of [^]: 
is: / is insensitive to [^] as judged by the corresponding value of aworthey (see text), mv: / 
is multi-valued, neg: I is negative (for those that are in linear Wx units). 

*I is a strong, well-behaved [^]-indicator over a broad [^] range down to XMP values (at 
least [^] < —5.0) for at least some GK star Teg values. 
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Table 3. LTE fitting function coefficients, Cn, for independent variables log 6* and [^] for 
stars of 3750 < < 6500 K at IDS spectral resolution for our nine Lick-XMP indices. 

Quantities in brackets below each fitted value are the 1-a “htting errors” (see text) 
estimated from and may be interpreted as 68% conhdence intervals. 


Polynomial term Lick index, / 


Term (n) 

Fe4383 

Fe4531 

Fe5015 

Mgi 

Mg b 

a 

5.7407 

5.1054 

8.5322 

0.0547 

0.7475 

log6' 

42.7841 

26.6568 

38.7751 

1.6149 

22.3066 

a 

(1.9465) 

(0.9536) 

(1.1653) 

(0.1791) 

(0.9649) 

log^ 6 

34.3084 

47.8340 

32.7787 

16.8150 

150.1230 

a 

(14.3478) 

(7.0287) 

(8.5895) 

(1.3200) 

(7.1125) 

log^e 

-383.6426 

6.9564 

-16.5491 

54.8095 

208.7209 

a 

(135.1847) 

(66.2243) 

(80.9304) 

(12.4368) 

(67.0145) 

[f] 

1.3678 

2.2633 

4.9326 

0.0345 

0.5357 

a 

(0.1858) 

(0.0910) 

(0.1112) 

(0.0171) 

(0.0921) 

w? 

-0.1555 

0.2680 

0.9472 

0.0021 

0.1526 

a 

(0.0749) 

(0.0367) 

(0.0448) 

(0.0069) 

(0.0371) 


-0.0376 

0.0050 

0.0600 

-0.0005 

0.0132 

a 

(0.0083) 

(0.0041) 

(0.0050) 

(0.0008) 

(0.0041) 


8.1137 

7.0252 

13.5876 

0.7520 

11.3125 

a 

(1.1698) 

(0.5731) 

(0.7003) 

(0.1076) 

(0.5799) 


0.2170 

0.3786 

1.1877 

0.0629 

1.2825 

a 

(0.1900) 

(0.0931) 

(0.1138) 

(0.0175) 

(0.0942) 


2.8481 

7.2659 

5.3031 

3.4235 

31.2653 

a 

(4.3128) 

(2.1128) 

(2.5820) 

(0.3968) 

(2.1380) 

Total a 

136.0316 

66.6392 

81.4374 

12.5147 

67.4343 

Red 

0.1755 

0.0421 

0.0629 

0.0015 

0.0431 

Term (n) 

Fe5270 

Fe5335 

Fe5406 

Na D 


C'o 

3.3993 

3.5991 

1.9507 

1.2374 


log 6* 

24.7325 

22.7161 

19.5786 

20.5914 


a 

(1.0332) 

(0.7873) 

(0.5449) 

(1.2010) 


\og^e 

43.8640 

27.8599 

55.6521 

151.3475 


a 

(7.6155) 

(5.8036) 

(4.0165) 

(8.8525) 


log^ 9 

-179.7178 

-169.1160 

-68.7232 

407.8087 
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Table 3—Continued 


Polynomial term Lick index, / 


a 

(71.7529) 

(54.6813) 

(37.8434) 

(83.4082) 

m 

1.7790 

1.8678 

1.1233 

1.0780 

a 

(0.0986) 

(0.0752) 

(0.0520) 

(0.1146) 


0.3000 

0.3264 

0.2176 

0.2884 

a 

(0.0398) 

(0.0303) 

(0.0210) 

(0.0462) 

[W 

0.0165 

0.0193 

0.0140 

0.0228 

a 

(0.0044) 

(0.0033) 

(0.0023) 

(0.0051) 

[f]log0 

6.3203 

6.8152 

6.9676 

12.0989 

a 

(0.6209) 

(0.4732) 

(0.3275) 

(0.7218) 


0.3882 

0.5507 

0.6483 

1.4020 

a 

(0.1009) 

(0.0769) 

(0.0532) 

(0.1172) 

[M] log2 0 

3.9294 

1.8886 

8.9946 

33.6358 

a 

(2.2892) 

(1.7445) 

(1.2073) 

(2.6610) 

Total a 

72.2025 

55.0239 

38.0804 

83.9307 

Red 

0.0494 

0.0287 

0.0137 

0.0668 
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Table 4. Same as Table [3l but for spectra of SDSS spectral resolution. 


Polynomial term Lick index, / 


Term (n) 

Fe4383 

Fe4531 

Fe5015 

Mgi 

Mg b 

C'o 

6.3207 

5.8713 

9.9981 

0.0596 

0.7159 

log^ 

46.5420 

29.7759 

47.4312 

1.6445 

21.1327 

a 

(2.1429) 

(1.1107) 

(1.0168) 

(0.1790) 

(1.0098) 

hg^e 

37.7593 

44.8405 

37.9836 

16.8348 

136.5346 

a 

(15.7958) 

(8.1873) 

(7.4950) 

(1.3194) 

(7.4430) 

hg^e 

-426.9500 

-45.0185 

-124.2286 

54.8555 

149.0733 

a 

(148.8280) 

(77.1404) 

(70.6182) 

(12.4315) 

(70.1279) 

1^1 

1.5700 

2.5478 

5.5798 

0.0376 

0.5076 

a 

(0.2046) 

(0.1060) 

(0.0971) 

(0.0171) 

(0.0964) 


-0.1475 

0.2888 

1.0275 

0.0028 

0.1489 

a 

(0.0825) 

(0.0427) 

(0.0391) 

(0.0069) 

(0.0389) 

[W 

-0.0392 

0.0041 

0.0624 

-0.0004 

0.0134 

a 

(0.0091) 

(0.0047) 

(0.0043) 

(0.0008) 

(0.0043) 

[f]log0 

8.5967 

7.4703 

15.1617 

0.7638 

10.7092 

a 

(1.2879) 

(0.6675) 

(0.6111) 

(0.1076) 

(0.6068) 


0.1880 

0.3720 

1.2213 

0.0641 

1.2324 

a 

(0.2092) 

(0.1084) 

(0.0993) 

(0.0175) 

(0.0986) 

[f]log^^ 

2.3983 

5.9506 

4.2027 

3.4272 

28.2502 

a 

(4.7481) 

(2.4610) 

(2.2530) 

(0.3966) 

(2.2373) 

Total a 

149.7603 

77.6237 

71.0606 

12.5094 

70.5673 

Red 

0.2127 

0.0571 

0.0479 

0.0015 

0.0472 

Term [n] 

Fe5270 

Fe5335 

Fe5406 

Na D 


Co 

3.6687 

4.2195 

2.2520 

1.3092 


log 9 

26.8280 

26.0111 

23.4383 

21.0015 


a 

(1.1715) 

(0.8517) 

(0.6402) 

(1.2074) 


log^ 9 

44.6660 

32.4921 

67.7726 

152.7714 


a 

(8.6353) 

(6.2780) 

(4.7191) 

(8.8997) 


log^9 

-205.3421 

-169.5044 

-80.2579 

412.8384 
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Table 4—Continued 


Polynomial term Lick index, / 


a 

(81.3620) 

(59.1509) 

(44.4633) 

(83.8526) 

m 

1.8747 

2.2036 

1.3076 

1.1254 

a 

(0.1118) 

(0.0813) 

(0.0611) 

(0.1152) 


0.3043 

0.3826 

0.2548 

0.2983 

a 

(0.0451) 

(0.0328) 

(0.0246) 

(0.0465) 

[W 

0.0158 

0.0222 

0.0164 

0.0235 

a 

(0.0050) 

(0.0036) 

(0.0027) 

(0.0051) 

[f]log^ 

6.5574 

8.0558 

8.4043 

12.2936 

a 

(0.7041) 

(0.5119) 

(0.3848) 

(0.7256) 


0.3685 

0.6646 

0.7870 

1.4220 

a 

(0.1144) 

(0.0831) 

(0.0625) 

(0.1179) 

[M] log2 0 

3.4790 

2.5911 

11.0017 

33.9557 

a 

(2.5957) 

(1.8871) 

(1.4185) 

(2.6752) 

Total a 

81.8717 

59.5215 

44.7419 

84.3780 

Red 

0.0636 

0.0336 

0.0190 

0.0675 
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Table 5. Same as Table [3l but for independent variables V — K and [^]. 


Polynomial term Lick index, / 


Term (n) 

Fe4383 

Fe4531 

Fe5015 

Mgi 

Mg b 

C^o 

2.0247 

2.0280 

0.0297 

0.4264 

1.9868 

V - K 

-2.8895 

-0.4509 

4.5572 

-0.6166 

-3.5266 

a 

(1.3722) 

(0.8606) 

(0.8227) 

(0.1854) 

(0.4996) 

V-K^ 

3.0594 

1.1931 

-0.3775 

0.2549 

1.5847 

a 

(0.5385) 

(0.3377) 

(0.3229) 

(0.0728) 

(0.1961) 

V -K^ 

-0.4494 

-0.1639 

0.0199 

-0.0236 

-0.1178 

a 

(0.0657) 

(0.0412) 

(0.0394) 

(0.0089) 

(0.0239) 

m 

-0.4489 

0.9816 

2.2984 

-0.0033 

-0.6624 

a 

(0.3320) 

(0.2082) 

(0.1991) 

(0.0449) 

(0.1209) 


-0.1629 

0.2169 

0.7254 

-0.0077 

-0.1158 

a 

(0.0792) 

(0.0497) 

(0.0475) 

(0.0107) 

(0.0288) 

[f\^ 

-0.0373 

0.0022 

0.0540 

-0.0012 

0.0042 

a 

(0.0077) 

(0.0048) 

(0.0046) 

(0.0010) 

(0.0028) 

mv-K 

0.9985 

0.6180 

1.2306 

-0.0144 

0.1333 

a 

(0.1784) 

(0.1119) 

(0.1069) 

(0.0241) 

(0.0649) 

[f\^V-K 

0.0059 

0.0097 

0.0696 

0.0009 

0.0749 

a 

(0.0158) 

(0.0099) 

(0.0095) 

(0.0021) 

(0.0057) 

SIS 

1 

to 

-0.0728 

-0.0308 

-0.0576 

0.0109 

0.1277 

a 

(0.0292) 

(0.0183) 

(0.0175) 

(0.0040) 

(0.0106) 

Total a 

1.5253 

0.9566 

0.9146 

0.2061 

0.5554 

Red 

0.0149 

0.0059 

0.0054 

0.0003 

0.0020 

Term (n) 

Fe5270 

Fe5335 

Fe5406 

Na D 


Co 

2.3521 

2.1055 

1.7732 

2.0387 


V-K 

-2.7555 

-2.0371 

-2.7317 

-2.8668 


a 

(0.8051) 

(0.5521) 

(0.5025) 

(0.9161) 


V-K^ 

2.0821 

1.7716 

1.7435 

1.2867 


a 

(0.3159) 

(0.2166) 

(0.1972) 

(0.3595) 


V -K^ 

-0.2818 

-0.2483 

-0.2162 

-0.0768 
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Table 5—Continued 


Polynomial term Lick index, / 


a 

(0.0386) 

(0.0264) 

(0.0241) 

(0.0439) 

m 

0.4714 

0.3801 

-0.0986 

-0.1966 

a 

(0.1948) 

(0.1336) 

(0.1216) 

(0.2217) 


0.2452 

0.2306 

0.0954 

0.0009 

a 

(0.0465) 

(0.0319) 

(0.0290) 

(0.0529) 

[|]^ 

0.0158 

0.0189 

0.0119 

0.0132 

a 

(0.0045) 

(0.0031) 

(0.0028) 

(0.0051) 

mv-K 

0.6728 

0.7581 

0.5120 

0.1566 

a 

(0.1046) 

(0.0718) 

(0.0653) 

(0.1191) 

W?V-K 

0.0212 

0.0406 

0.0440 

0.0805 

a 

(0.0092) 

(0.0063) 

(0.0058) 

(0.0105) 

[M]!^ - 

-0.0390 

-0.0413 

0.0054 

0.1307 

a 

(0.0172) 

(0.0118) 

(0.0107) 

(0.0195) 

Total a 

0.8949 

0.6137 

0.5586 

1.0183 

Red 

0.0051 

0.0024 

0.0020 

0.0067 
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Table 6. LTE Partial I derivatives with respect to and [^], at each grid value of Tes 
and [^], in index units and in units of a for SDSS spectral resolution for our nine 
Lick-XMP indices. Sample entries only - the entire table is available electronically. 


Lick index, / 


TeS 


Index 

[f] 

6500 

6000 

5500 

5250 

5000 

4750 

4500 

4250 

Fe4383 

-0.0 

-0.1508 

-0.2423 

-0.3303 

-0.3706 

-0.4064 

-0.4355 

-0.4551 

-0.4610 

Fe4383 

-0.5 

-0.1241 

-0.2129 

-0.2975 

-0.3358 

-0.3694 

-0.3961 

-0.4129 

-0.4157 

Fe4383 

-1.0 

-0.0981 

-0.1841 

-0.2654 

-0.3017 

-0.3332 

-0.3575 

-0.3717 

-0.3714 

Fe4383 

-1.5 

-0.0727 

-0.1560 

-0.2340 

-0.2685 

-0.2978 

-0.3198 

-0.3313 

-0.3281 

Fe4383 

-2.0 

-0.0480 

-0.1286 

-0.2034 

-0.2360 

-0.2633 

-0.2830 

-0.2919 

-0.2857 


Lick index, I 



100Kx,g^ 

|[M/H]/0- 





Index 

[|] 

6500 

6000 

5500 

5250 

5000 

4750 

4500 

4250 

Fe4383 

-0.0 

-0.0010 

-0.0016 

-0.0022 

-0.0025 

-0.0027 

-0.0029 

-0.0030 

-0.0031 

Fe4383 

-0.5 

-0.0008 

-0.0014 

-0.0020 

-0.0022 

-0.0025 

-0.0026 

-0.0028 

-0.0028 

Fe4383 

-1.0 

-0.0007 

-0.0012 

-0.0018 

-0.0020 

-0.0022 

-0.0024 

-0.0025 

-0.0025 

Fe4383 

-1.5 

-0.0005 

-0.0010 

-0.0016 

-0.0018 

-0.0020 

-0.0021 

-0.0022 

-0.0022 

Fe4383 

-2.0 

-0.0003 

-0.0009 

-0.0014 

-0.0016 

-0.0018 

-0.0019 

-0.0019 

-0.0019 


Lick index, I 




0 5 X 

U.O X 

I'Lff 





Index 

[|] 

6500 

6000 

5500 

5250 

5000 

4750 

4500 

4250 

Fe4383 

-0.0 

0.3247 

0.4664 

0.6237 

0.7092 

0.7999 

0.8964 

0.9994 

1.1098 

Fe4383 

-0.5 

0.3942 

0.5326 

0.6863 

0.7699 

0.8586 

0.9530 

1.0539 

1.1619 

Fe4383 

-1.0 

0.4342 

0.5693 

0.7195 

0.8012 

0.8879 

0.9803 

1.0789 

1.1846 

Fe4383 

-1.5 

0.4449 

0.5767 

0.7233 

0.8031 

0.8879 

0.9781 

1.0745 

1.1779 

Fe4383 

-2.0 

0.4261 

0.5547 

0.6977 

0.7756 

0.8584 

0.9465 

1.0408 

1.1418 
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Table 6—Continued 


Lick index, / 


TeS 


Lick index, I 


0.5 X ^^\ng/a 


Index 

[f] 

6500 

6000 

5500 

5250 

5000 

4750 

4500 

4250 

Fe4383 

-0.0 

0.0022 

0.0031 

0.0042 

0.0047 

0.0053 

0.0060 

0.0067 

0.0074 

Fe4383 

-0.5 

0.0026 

0.0036 

0.0046 

0.0051 

0.0057 

0.0064 

0.0070 

0.0078 

Fe4383 

-1.0 

0.0029 

0.0038 

0.0048 

0.0053 

0.0059 

0.0065 

0.0072 

0.0079 

Fe4383 

-1.5 

0.0030 

0.0039 

0.0048 

0.0054 

0.0059 

0.0065 

0.0072 

0.0079 

Fe4383 

-2.0 

0.0028 

0.0037 

0.0047 

0.0052 

0.0057 

0.0063 

0.0069 

0.0076 
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Table 7. Same as Tabled but with respect to V — K and [^], at each grid value oiV — K 
for [^] = 0.0 models and [^]. Note that V — K is [^]-dependent and the V — K values 
used to label the columns are only valid for [^] = 0.0 (Fig. H] gives an indication of the 
[^(-dependence on 1/ — K.) Sample entries only - the entire table is available electronically. 


0.25 mag x q(^v-k) \wh] 
Lick index, / {V - K) at [M/H] = 0.0 


Index 

[f] 

0.887 

1.211 

1.578 

1.792 

2.017 

2.290 

2.585 

Fe4383 

-0.0 

0.6126 

0.8463 

1.0189 

1.0744 

1.0969 

1.0748 

0.9900 

Fe4383 

-0.5 

0.5123 

0.7459 

0.9227 

0.9801 

1.0083 

0.9960 

0.9239 

Fe4383 

-1.0 

0.4144 

0.6498 

0.8299 

0.8884 

0.9203 

0.9136 

0.8495 

Fe4383 

-1.5 

0.3137 

0.5543 

0.7361 

0.7977 

0.8326 

0.8309 

0.7770 

Fe4383 

-2.0 

0.2136 

0.4555 

0.6419 

0.7063 

0.7449 

0.7487 

0.7029 


Lick index, I 


0.25 mag 


Index 

[f] 

0.887 

1.211 

1.578 

1.792 

2.017 

2.290 

2.585 

Fe4383 

-0.0 

-2.8131 

-2.2252 

-1.6230 

-1.3032 

-0.9918 

-0.6481 

-0.3189 

Fe4383 

-0.5 

-2.9024 

-2.3225 

-1.7227 

-1.4147 

-1.1006 

-0.7738 

-0.4398 

Fe4383 

-1.0 

-2.9874 

-2.4089 

-1.8095 

-1.5086 

-1.1932 

-0.8744 

-0.5388 

Fe4383 

-1.5 

-3.0797 

-2.4949 

-1.9022 

-1.5992 

-1.2884 

-0.9726 

-0.6457 

Fe4383 

-2.0 

-3.1715 

-2.5911 

-1.9978 

-1.6955 

-1.3860 

-1.0757 

-0.7484 


Lick index, I 



0.5 X gjMjl 

iV-K) 




Index 

If) 

0.887 

1.211 

1.578 

1.792 

2.017 

2.290 

2.585 

Fe4383 

-0.0 

-1.0540 

-0.8495 

-0.6315 

-0.5151 

-0.3944 

-0.2651 

-0.1258 

Fe4383 

-0.5 

-0.9953 

-0.7911 

-0.5733 

-0.4570 

-0.3365 

-0.2074 

-0.0683 

Fe4383 

-1.0 

-0.9658 

-0.7618 

-0.5442 

-0.4281 

-0.3077 

-0.1787 

-0.0398 

Fe4383 

-1.5 

-0.9654 

-0.7616 

-0.5443 

-0.4283 

-0.3080 

-0.1792 

-0.0405 

Fe4383 

-2.0 

-0.9940 

-0.7905 

-0.5734 

-0.4575 

-0.3374 

-0.2088 

-0.0702 
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Table 7—Continued 


0.25 mag x Q^v-K) \wm 
Lick index, I {V - K) at [M/H] = 0.0 


Lick index, / 


0.5 X 


Index 

[f] 

0.887 

1.211 

1.578 

1.792 

2.017 

2.290 

2.585 

Fe4383 

-0.0 

-0.7299 

-0.5883 

-0.4373 

-0.3567 

-0.2731 

-0.1836 

-0.0871 

Fe4383 

-0.5 

-0.6893 

-0.5479 

-0.3970 

-0.3165 

-0.2330 

-0.1436 

-0.0473 

Fe4383 

-1.0 

-0.6689 

-0.5276 

-0.3769 

-0.2965 

-0.2131 

-0.1238 

-0.0276 

Fe4383 

-1.5 

-0.6686 

-0.5274 

-0.3769 

-0.2966 

-0.2133 

-0.1241 

-0.0280 

Fe4383 

-2.0 

-0.6884 

-0.5474 

-0.3971 

-0.3169 

-0.2337 

-0.1446 

-0.0486 








